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INTRODUCTION 

The improving precision of experiments in the high energy physics motivates the- 
oretical studies of quantum corrections to various processes. In the two-loop approx- 
imation this is connected with great computational difficulties, especially if there are 
several mass scales involved in the process, which is a typical situation in the case 
of electroweak or mixed chromodynamic and electroweak corrections. Recently a new 
method has been proposed for the evaluation of scalar two-loop vertex and propagator 
functions |TJ, . It has also been shown that a similar approach works even for the four- 
point functions [Q. In this talk I present a few examples which illustrate the principle 
of this method. 

The aim will be to obtain a double integral representation which is suitable for 
numerical evaluation. In the following section I will derive it for a special case of the 
vertex function with zero momentum transfer. The next section refers to the general 
case of a planar vertex function with space-like values of external momenta, and the last 
one shows an example of dealing with ultraviolet divergent diagrams. The examples of 
two- loop functions to be considered in this paper are depicted in Fig. l(a,b,c). 



VERTEX FUNCTION AT ZERO MOMENTUM TRANSFER 

While ref. |§ describes the general method of computing the two-loop vertex func- 
tion, here the method is illustrated with the special case of zero-momentum transfer 

1 Submitted to the proceedings of the NATO Advanced Study Institute "Frontiers in Particle 
Physics" , Cargese, August 1994. 
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putation becomes mucn simpler ana it is easy to write aown explicit lormuias. ine 
diagram and numbering of lines is depicted in Fig. 1(b). The four momenta in the rest 
frame of the external particle are (p, are outgoing) 



P% = -rf = (q, 0,0,0), 
F = (l+ + l-,l+-l_X 
P = (k + + k_,k + -k-,k±), (1) 

and the momentum k runs along the lines (4,3,5,6) and / along (1,3,2). We shall see 
later that one obtains different, but equivalent, formulas if one uses a different routing 
of the momenta. The two- loop function of Fig. 1(b) is 



With the present choice of momenta, and with s = l\ and t = kj_, the explicit form of 
the propagators is 

Pi,2 = (2l + -q)(2L-q)-s-ml >2 +ir), 

P 3 = 4(i + + k + )(l_ + k_)-s-t-2Vsiz-m 2 3 + iri, 

P 4 , 5 = {2k+ + q)(2k- + q)-t-ml & + iTi, 

P 6 = 4k + k_ -t-ml+ir). (3) 

There is only one propagator, P 3 , through which both internal momenta flow, and z 
denotes the cosine of the angle between the two perpendicular momentum vectors k± 
and l± . Integrations over the two angular variables describing the absolute and relative 
configuration of k± and l± can now be performed and we obtain from (§) 

r 11 
V (q 2 ) = -An 2 / dk+dl+dsdtdk-dl- , , (4) 

with 

A = 4(1+ + k+)(l- + k-) - s -t-ml + ir], 
B 2 = 4st. (5) 

The integrations over the and /_ are done with help of contour integrals. It turns 
out that the singularities in P4 and P5 do not contribute, while Pi, P2 and Pq contribute 
only for k + and l + lying in a triangular region T in the k+l+ plane delimited by the 
lines 

k + + l + = 0, 
k+ = 0, 

h - |. (6) 

The function Vn becomes now 



Jo 
x V 1 1 



(7) 
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_ t + mj _ s + ml 2 q 

h ~ lifeT' h > 2 ~ 2(2Z + - g) + 2" (8) 

Substituting the explicit formulas for the propagators one obtains 



X 



with 



q 2 (m\ - ml) J Jt + + 2l + -qJo (t + t 4 )(t + t 5 ) 

Hdsl - 1 - 1 =), (9) 

J ° \ ^J(at + b 2 + cs) 2 - Ast yj(at + h + cs) 2 - Ast ) 



it 



2fc + + g 
2Z + -g' 

k 5 = f t {( 2 ^+^)(^+^)-< 5 }- ( 10 ) 

For the sake of simplicity let us assume that q 2 lies below all thresholds so that the 
function V is real. In such case the integration over s is elementary 

1 1 




(at + b 2 + cs) 2 - Ast J (at + bi + cs) 2 - Ast 



l_ ln t(l-ca)-cb 2 

c t(l - ca) -cbx y ' 



and in the integration over t one encounters dilogarithms, 

ln|l -y 



Li 2 (x 



y 

The final result is 

v( 2 )= 4tt 4 ro dk+ ri/ 2 

q(vn\ — m 2 )(m 2 — m\) J-q/2 2k + + q J-k+ 



(13) 



with ti j2 = —061,2/(1 — ca). 

The complication which arises in the case of the non-zero momentum transfer con- 
sists in the fact that the propagators after substitution of the appropriate value 
for do not have the simple form ^(m 2 . — m\) but retain dependence on the variables 
l + and s. This leads to a more complicated form for the integrations over s and t, and 
the final formula contains dilogarithms as well as Clausen functions. There are also 
two additional residues which contribute and each contribution in general comes from 
a different triangle in the k + l + plane. 
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to nave an alternative tormuia wmcn provides a cross cneck ana a test 01 accuracy. 
Such formula can be derived by choosing the internal momenta in such way that k runs 
through the lines (1, 2, 5, 6, 4) and I through (1, 2, 3). 
With this choice the propagators are 



Pl,2 
Pi 

-P-4,5 
P 6 



(2k+ + 21+ + q){2k_ + 2Z_ + q) - s 
4Z + L_ — s — 772,3 + 
(2k + + q)(2k^ + q) - t - m 
4fc + fc_ — t — m\ + irj. 



t - 2\[stz 



rn 



1)2 



4,5 



+ irj, 



(14) 



There are now two propagators which depend on z: P\ and P%, but only one combination 
of propagators (P 3 and P 6 ) whose singularities contribute to the contour integrations 
over and /_ . The triangular region of the integration over k + and l + is now limited 
by the lines 



Q 
2 



0. 
0. 



(15) 



After the integration over the angular variables and over k_ and /_ we get 




(16) 



*k 6 ,l-->l 3 



with 



t + m 




m 



1,2- 



(17) 



The integrations over s and t proceed in exactly the same way as in the previous 
calculation. Finally, we can make the shift l + — > l + — q/2. It turns out that this change 
of variables not only makes the region of integration in the k + l + plane equal to the 
triangle T defined in @, but the whole formula for Vi becomes almost the same as 
formula (pf), the only difference being the coefficients bi, which in the present case are 



b'=2(k + + l + ) 



777« 



2k, 



+ 



2L 



q 



+ q — m 



1,2- 



The equivalence of the two formulas can be checked after integrating over k + and l + . 
It provides an excellent cross check for the numerical calculation. 

In practical calculations one can encounter a mass configuration in which m\ = m 2 . 
In this case the formula simplifies: 



V (q 2 ,m 1 = m 2 ) 



4tt 4 



x 



dk A 



q/2 



q/2 
-h+ 



q 2 (ml — m 
h 

dh 



k+ + L 



1 . h 1 

In 

£5 — ^0 ^4 — 



In 



u 



(19) 
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SPACE-LIKE EXTERNAL MOMENTA 



If the external momenta have space-like values the computation of the propagator 
and vertex diagrams is greatly simplified since the internal particles do not become 
on-shell. In particular we can easily check the analytical results obtained for the two- 
point and the planar three-point functions when all internal particles are massless. 
While the result for the two-point function (see Fig. 1(a)) has been know for long 
time H], the much more complex formulas for the vertex functions (of both planar 
and crossed topologies) have been obtained only very recently 0. We present here 
numerical evaluation of the vertex function with all internal masses equal m and space- 
like external momenta. In the limit m^Owe reproduce the result of ||. 

For the numerical calculation it is convenient to choose such reference frame that 
the external outgoing momenta become (according to the notation of Fig. 1(b)) 

Pi = (e, gi,0,0), 
P2 = (-e,g 2 ,0,0), 

P£ = (0,-gi-g 2 ,0,0). (20) 

Repeating the calculations described in the previous section we arrive at a double 
integral representation which is easy to evalulate numerically. Fig. 2 shows the ratio of 
the vertex function 

v<*.AM = Jf n££w (21) 

to the value of the vertex at zero internal masses 

^A P ^(f)V>(Jj), (22) 



where the function $( 2 ) has been derived in ||. In formula ( |2~T| ) Pj denote propagators 
defined analogously to the formula ©). For the purpose of numerical calculation we 
choose one arbitrary configuration of external momenta p\ = —1, p\ = —4, p\ = —25. 
We see that for very small masses the ratio of the two formulas becomes unity which 
confirms the analytical result of Ussyukina and Davydychev. 



DIVERGENT INTEGRALS 

The method of calculation of two-loop diagrams described here is limited to the 
four-dimensional space. In dealing with divergent integrals we first have to find another 
diagram with the same divergent part but simple enough to be computed analytically. 
The difference of the two diagrams can then be calculated in four dimensions and the 
final result is obtained by adding the analytical formula for the simpler diagram. Such 
procedure, based on the representation of the two-point functions proposed in J7|, has 
been described in ||. In the present section I illustrate an analogous procedure in the 
framework of the representation which works for both two- and three-point functions, 
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propagators {rig. i(cjj. 

Considerable effort has been recently devoted to the investigation of this diagram. 
Asymptotic expansions have been derived in the papers || [T(|, and explicit expres- 
sions in terms of generalized hypergeometric, or Lauricella, functions were obtained in 
ref . [|TTJ . The same diagram has also been analyzed in jl2] . It has been noted in [ 13| that 
a general two-loop diagram can formally be expressed as a sunrise diagram with masses 
and the external momentum being functions of Feynman parameters over which one 
can integrate numerically The latter reference gave a very convenient one-dimensional 
integral representation for this diagram. 

The value of the sunrise diagram is 



-/E\2ul 



S(p 2 ,m 1 ,m 2 ,m 3 ) = (tt< 

with 7s being Euler's constant and 

Pi 
P2 

P 3 = k 2 -ml + 



II dDkdD 'w 



P2P3 



(23) 



(I + k + p) 2 — m\ + ir], 



I 2 — m\ + if], 



17], 



(24) 



and since the sunrise diagram is ultraviolet divergent we have to compute it in D = 
4 — 2uj dimensions. 

It has been shown in the previous sections that the triangular regions over which one 
has to perform the final two integrations numerically are determined only by the values 
of external momenta, and are independent of masses of particles inside the diagram. 
Therefore it is convenient to choose for the subtraction a diagram which differs from the 
diagram we are interested in only by the values of internal masses. In the present case 
we choose a diagram with vanishing m 2 and m 3 which can be computed analytically 

S(p 2 ,m u 0,0) ' 



— IT 



— ? + — (1 -21nm? 



+ In ml (in 



7T 
1 

-,2 




oh 

(25) 



and the value of a diagram with arbitrary masses can be expressed by 

S(p 2 ,m 1 ,m 2 ,m 3 ) = S(p 2 , m Xy 0, 0) + A(p 2 ,m 1 ,m 2 ,m 3 ) (26) 

where A(p 2 , mi, m 2 , m 3 ) = A is free from divergences and can be computed using our 
method. For simplicity we only consider the case of p 2 < m\ where both diagrams are 
real. 

After the integration over angular variables as in (d) and over fc_ and with help 
of contour integrals we obtain 



a - a 1 / / d ^4±[r d sdt 

T k + l + J Jo 



A 



(^42 _ 52)3/2 (Ai-B 2 yi 2 i 



(27) 



with 



.4 



(2l + + 2k ++P ) (^ + ^+P 



2L 



2k, 



— m l — s — t + irj 



at + b + cs 



(28) 
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integration is a triangle 1 delimited by tne lines fc + = U, < + = U and l + + k + = —p/Z. 
The integrations over s and t are easy 



I 



00 j at + b + cs 1 

ds 5^ = — for a, 6, c < (29) 

[(at + 6 + csf - Astf 2 (l-ac)t-bc 
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and finally we arrive at 



p/2 



k+ p(2l+ + 2k + + p) ( p+ |^ + |£) { 2i + + 2k + +p)-ml 

(30) 



Thus we have found a double integral representation of the sunrise diagram. One of 
the k + , l + integrations can still be carried out, and since the argument of the logarithm 
is a polynomial of the second degree the result will in general involve dilogarithms of 
complex arguments even below the threshold. For the purpose of numerical evaluation 
it may be convenient to work with a double-integral, but explicitly real representation. 



ACKNOWLEDGMENTS 

I thank D. Broadhurst, K.G. Chetyrkin, and A.I. Davydychev for discussion and 
advice, and B. Krause and M. Steinhauser for checking some of the formulas. I am 
very grateful to the organizers of the Cargese Summer Institute for the opportunity 
to take part in this great event. I thank Graduiertenkolleg Elementarteilchenphysik at 
the University of Karlsruhe for support. 

References 

D. Kreimer, Phys. Lett. B292 (1992) 341. 



[2 

[3] 
[4] 
[5] 
[6] 

[7 

[8 

[9 
[10; 



A. Czarnecki, U. Kilian, and D. Kreimer, New representation of two-loop propaga- 
tor and vertex functions, |hep-ph / 9405423| , in press in Nucl. Phys. B. 



D. Kreimer, A short note on two-loop box functions, |iep-ph/9407234 
J.L. Rosner, Ann. Phys. 44 (1967) 11. 

K.G. Chetyrkin and F.V. Tkachov, Nucl. Phys. B192 (1981) 159. 

N.I. Ussyukina and A.I. Davydychev, Phys. Lett. B298 (1993) 363; Yad. Fiz. 56 
(1993) 172; Phys. Lett. B332 (1994) 159. 

D. Kreimer, Phys. Lett. B273 (1992) 277. 

FA. Berends and J.B. Tausk, Nulc. Phys. B421 (1994) 456. 

A.I. Davydychev and J.B. Tausk, Nucl. Phys. B397 (1993) 123. 

A.I. Davydychev, VA. Smirnov, and J.B. Tausk, Nucl. Phys. B410 (1993) 325. 



7 



b. oauDerger, t .A. aerenas, M. Donm, ana M. tfuza, |tiep-pn/^4U^^aq . 

[12] F.A. Lunev, On evaluation of two-loop self-energy diagram with three propagators, 
|hep-th/9408T6Tj 

[13] A. Ghinculov and J.J. van der Bij, Massive two-loop diagrams: the Higgs propa- 
gator, [hep-ph/ 940541 8 . 



8 



